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"Hot Jupiter" extrasolar planets are expected to be tidally locked because they are close 
(<0.05 astronomical units, where 1 AU is the average Sun-Earth distance) to their parent 
stars, resulting in permanent daysides and nightsides. By observing systems where the planet 
and star periodically eclipse each other, several groups have been able to estimate the tem- 
peratures of the daysides of these planets 1-3 . A key question is whether the atmosphere is 
able to transport the energy incident upon the dayside to the nightside, which will deter- 
mine the temperature at different points on the planet's surface. Here we report observa- 
tions of HD 189733, the closest of these eclipsing planetary systems 4-6 , over half an orbital 
period, from which we can construct a 'map' of the distribution of temperatures. We de- 
tected the increase in brightness as the dayside of the planet rotated into view. We estimate 
a minimum brightness temperature of 973 ±33 K and a maximum brightness temperature of 
1212±11 K at a wavelength of 8 /xm, indicating that energy from the irradiated dayside is 
efficiently redistributed throughout the atmosphere, in contrast to a recent claim for another 
hot Jupiter 7 . Our data indicate that the peak hemisphere-integrated brightness occurs 16±6 
degrees before opposition, corresponding to a hot spot shifted east of the substellar point. 
The secondary eclipse (when the planet moves behind the star) occurs 120±24 s later than 
predicted, which may indicate a slightly eccentric orbit. 

We monitored HD 189733 continuously over a 33.1 hour period using the 8 /im channel of 
the InfraRed Array Camera (IRAQ 8 on the Spitzer Space Telescope 9 , observing in subarray mode 



1 



with a cadence of 0.4 s. Our observations spanned slightly more than half of the planet's orbit, 
beginning 2.6 hours before the start of the transit (when the planet moves in front of the star) and 
ending 1.9 hours after the end of the secondary eclipse. This gave us a total of 278,528 32 x 32 
pixel images. We found that there was a gradual detector-induced rise of up to 10% in the signal 
measured in individual pixels over time. This rise is illumination-dependent; pixels with high 
levels of illumination (greater than 250 MJy sr" 1 ) converge to a constant value within the first two 
hours of observations and lower-flux pixels increase linearly over time. We characterize this effect 
by producing a time series of the signal in a series of annuli of increasing radius centered on the star 
(masking out a 5-pixel- wide box centered on HD 189733's smaller, fainter M dwarf companion 10 ). 
This set of curves describes the behavior of the ramp for different illumination levels. 

To correct our images, we estimate the median illumination for each pixel in the array, and 
interpolate over our base set of curves (scaling as the natural log of the illumination) to calculate a 
curve describing the behavior of that pixel. We correct for this instrumental effect by dividing the 
flux in each pixel in a given image by the value of the interpolated curve. Pixels with illumination 
levels higher than 210 MJy sr -1 are not corrected, as these pixels converge to a constant value 
before the transit. We subsequently measure the flux from the M dwarf companion and find that 
it is constant at a level of < 0.05%, indicating that the detector effect has been removed. We 
then use aperture photometry with a radius of 3.5 pixels to create our time series (see Fig. 1 for 
additional details). We chose the smallest aperture possible while still avoiding flux losses from 
the shifting position of the star on the array; only 33% of the total flux in our aperture comes from 
corrected pixels. We test the effect of our correction on the size of the observed signal by inserting 
an artificial phase variation signal into the images before applying our corrections; we find that the 
amplitude of this signal is reduced by only 13% of its total size. 

In our final time series (Fig. 1), we see a distinct rise in flux beginning shortly after the end 
of the transit and continuing until a time just prior to the beginning of the secondary eclipse. (The 
transit and secondary eclipse occur at orbital phases ~ and ~ 0.5, respectively.) We estimate its 
amplitude by fitting a small region of the phase curve around the peak with a quadratic function, 
and taking the maximum value of that function as the peak of the phase curve. After similarly fitting 
the region around the minimum, we estimate the total amplitude of this rise to be 0.12 ± 0.02%, 
with uncertainties that are dominated by our correction for the detector effect (we propagate this 
systematic effect in all of our stated uncertainties below). By comparing this variation to the depth 
of the secondary eclipse, we find that the minimum hemisphere-integrated flux is 64 ± 7% of the 
maximum flux. The peak in flux occurs 2.3 ± 0.8 hours before the center of the secondary eclipse, 
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corresponding to a position 16 ± 6 degrees before opposition. A possible confusing effect results 
from the fact that HD 189733 is an active star, with spots that cause the flux to vary by as much as 
±1.5% in visible light over its 13.4 day rotation period 6 . We estimate the size of this effect at 8 yum 
by treating both the star and the spots as blackbodies with temperatures of 5050 K and 4050 K 
respectively and scaling the variations observed at visible wavelengths to the appropriate 8 /im 
amplitude. Projecting these variations forward in time, we find that there could be a linear increase 
in flux of 0.1% over the period of our observations as the spots rotate into view. Importantly, we 
note that accounting for these spots would serve only to reduce the amplitude of the observed phase 
variation. As the shape of the observed variation is consistent with a genuine variation in the flux 
from the planet, we treat it as such in the discussions below. 

The high quality of our data allows us to derive more precise estimates of the parameters 
for the planetary system than are currently available 4 6 (see Fig. 2 and its legend for details). We 
calculate a brightness temperature of 1205. 1 ±9.3 K for the dayside of the planet from the observed 
depth of the secondary eclipse. We estimate that the planet has a minimum hemisphere-averaged 
brightness temperature of 973 ± 33 K occurring 6.7 ± 0.4 hours after the transit, and a maximum 
hemisphere- averaged brightness temperature of 1212 ± 11 K occurring 2.3 ± 0.8 hours before the 
onset of the secondary eclipse. 

We find the center of the transit occurs at tj = 2454037.611956 ± 0.000067 HJD (6 s error), 
while the center of the secondary eclipse occurs at time t u = 2454038.72294 ± 0.00027 HJD 
(24 s error), where the errors have been estimated using both a bootstrap method and a 10 5 step 
Markov chain. These are the most precise timing measurements of a transit and secondary eclipse 
to date. The transit occurs at the predicted time 6 , but the secondary eclipse occurs 150 ± 24 s 
later than its predicted time of half an orbital period after the transit. Because we observe both 
eclipses and the period is well-constrained, we are able to predict the time of secondary eclipse 
with no significant uncertainty. Part of the delay of the secondary eclipse is due to the light travel 
time across the system 13 of 30 s. The remaining delay is possibly due to (a) non-uniformity in the 
planet emission 1415 ; (b) third bodies in the system; or (c) an eccentric orbit. 

To estimate the magnitude of the first effect, we fit the observed phase variation with a simple 
model of the planet consisting of twelve longitudinal slices of constant brightness. The resulting 
model light curve is shown in Fig. lb, and the best-fit longitudinal flux values are shown in Fig. 3b 
(see figure legend for more details). Figure 3b is effectively a coarse 8 /im map of the planet with a 
resolution of 30 degrees in longitude and no resolution in latitude. Figure 3a shows this brightness 
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distribution projected onto the surface of the planet with an additional sinusoidal dependence on 
latitude included. Because we observe the planet over only half an orbital period, the error bars 
are largest for longitudes near 90 degrees west of the substellar point. Although this brightness 
distribution is a good fit for the later part of the phase curve (Fig. 1), a deviation is apparent 
near the transit; this fit could be improved by using a finer longitude resolution. We find that 
the brightest slice on the planet is 30 degrees east of the sub-stellar point. The faintest slice of 
the planet also (surprisingly) appears in the eastern hemisphere, 30 degrees west of the antistellar 
point. The brightest slice of the planet is roughly twice as bright as the faintest slice, corresponding 
to a temperature difference of ~350 K. This non-uniform brightness distribution changes the shape 
of the ingress and egress 1415 . Treating the planet as a uniformly bright disk in our fit creates an 
artificial delay of at most 20 s in the time of secondary eclipse. Thus, the planet's non-uniform 
emission cannot account for the 120 s delay of the secondary eclipse. 

This offset is unlikely to be the result of perturbations to the planet's mean motion from a 
third body in the system; such perturbations would shift the time of the transit as well, and we see 
no evidence for such a shift. This leaves the third option as the most likely explanation. If the time 
delay is attributed to eccentricity e, then ecosw = 0.0010 ± 0.0002 where w is the longitude of 
pericentre, indicating that the eccentricity is extremely small, but non-zero. This is surprising, as 
the time scale for orbital circularization is significantly shorter than the ages of these systems 1617 . 
This eccentricity is too small to have been detected by radial velocity measurements 4 6 . The 
observed delay is moderately inconsistent with the timing of the 16 /im eclipse 3 , which occurs 
29 ± 65 s later than predicted 6 . 

Atmosphere models allow us some insight into the factors that control the day-night temper- 
ature contrast. The response of a planet to stellar irradiation depends on a comparison between the 
radiative time scale (over which starlight absorption and infrared emission alter the temperature) 
and the advection time scale (over which air parcels travel between day and night sides) 18 20 . If 
the radiative time is much shorter than the advection time, the hot day side reradiates the absorbed 
stellar flux and the night side remains cold. If the radiative time greatly exceeds the advection time, 
however, then efficient thermal homogenization occurs. Radiative transfer models of highly irra- 
diated giant planets 19 23 predict that the bulk of absorption of stellar flux and emission of thermal 
flux occurs at pressures from tens of millibars to several bars, where the predicted radiative time 
scales 20 range from 10 4 -10 5 s. Advection times are less well constrained, but estimates of wind 
speeds 18 ' 24-28 (hundreds to thousands of m s _1 ) suggest advection times of ~ 10 5 s. Thus, current 
models suggest that the radiative time scale is comparable to the advective time scale, and temper- 
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ature differences could reach 1000 K. In contrast, the small flux variation that we observe implies 
that the time scale for altering the temperature by radiation modestly exceeds the time scale for 
homogenizing the temperature between the day and night sides. 

It is possible that the observed planetary flux emerges from deeper in the atmosphere than 
expected, where the radiative time scales are longer. In the 8 /im band, models suggest that H 2 
dominates the opacity, with additional contributions from CH 4 and collision-induced absorption of 
H 2 . Silicate cloud opacity is not expected at these temperatures 29 . If the radiative time constants 
are as small as expected 20 then supersonic wind speeds exceeding ~10 km s _1 (~ 4 times the 
sound speed) would be necessary to transport energy to the night side. The times of minimum 
and maximum flux also provide information on the planet's meteorology. Our observation that 
the minimum and maximum do not occur at phases of and 0.5, respectively, indicates advection 
of the temperature pattern by atmospheric winds 1 8 ' 24 ~ 28 ' 30 . The existence of a flux minimum and 
maximum on a single hemisphere suggests a complex pattern not yet captured in current circulation 
models. 

In contrast to the 8 /im phase variation for HD 189733b presented here, the 24 /im variation 
reported 7 for the non-transiting planet v And b was quite large. The reasons for the differing results 
are not immediately clear, although the sparse data sampling and unknown radius for v And b mean 
that the uncertainty in the inferred day-night contrast is much larger. A higher opacity at 24 /im and 
a lower surface gravity for v And b could lead to a photospheric pressure two times smaller, but this 
difference is likely insufficient to explain the discrepancy. The day side of v And b receives 50% 
more flux from its star, but it is unclear how this would affect the day/night temperature contrast. 
Secondary eclipse depths for several planets have been in good agreement with the predictions 
from simple ID models 19 ' 2123 that assume a uniform day-night temperature, consistent with our 
conclusions for HD 189733b. Taken together, these results argue for atmospheres that are very 
dark at visible wavelengths, likely absorbing 90% or more of the incident stellar flux, and at the 
same time able to transport much of this energy to the night side. 
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Figure 1: Observed phase variation for HD 189733b, with transit and secondary eclipse visible. 
We determine the location of the star by taking the weighted average of the flux contained in a 7- 
pixel box centered on the peak of the point spread function. We find that the shape of the observed 
variation is consistent for apertures between 3.5 — 7 pixels. IRAC takes images in sets of 64, and 
we found that the average fluxes in the first three images and the 58th image were consistently 
low. Additionally, 2% of the images had corrupted pixels within our aperture. We chose to trim 
both sets of images from our final time series. We also exclude the first 1.8 hours of data from 
our analysis, as our correction was designed to correct the data only beyond the start of the transit. 
We estimate the background flux by fitting a Gaussian function to a histogram of the fluxes from a 
subset of pixels located in the corners of the image. This background contributes 1.3% of the total 
flux in our aperture, and we subtract a constant value from our time series. The scatter in the final 
timeseries is 20% higher than predicted from photon noise alone; we use the standard deviation 
of the points after the end of the secondary eclipse as our error for each point. The stellar flux as 
measured at the center of the secondary eclipse is normalized to unity (dashed line), and the data 
is binned every 500 points (200 s). a and b show the same data, but in b the y axis is expanded to 
show the scale of the variation. The solid line is the phase curve for the best-fit model (Fig. 3). 

Figure 2: Time series of the transit and secondary eclipse. Data is binned every 100 points (40 s), 
with the out-of-transit fluxes normalized to unity, a shows the transit and b shows the secondary 
eclipse with best-fit eclipse curves overplotted, including timing offset, and residuals for the transit 
(c) and secondary eclipse (d) are plotted below. The out of transit data for the eclipses are nor- 
malized using a constant; we find the transit occurs 1 s earlier and the secondary eclipse occurs 
8 s earlier if we use a linear fit instead, an insignificant difference. We fit both eclipses, fixing the 
mass of the star 4 and allowing the transit times to vary freely 11 . From the primary eclipse, we find 
the radius of the planet is 1.137 ± 0.006 Rj up , the orbital inclination is 85.61 ± 0.04 degrees, and 
the radius of the star is 0.757 ± 0.003 Rs un ', the planet-star radius ratio is 0.1545 ± 0.0002. The 
formal uncertainty in the mass of the star introduces an additional error of ±1.8% in our estimates 
for the two radii. The depth of the secondary eclipse is 0.3381 ± 0.0055% in relative flux. Using 
a model 12 we predict a stellar brightness temperature of 4512 K in the 8 /im Spitzer bandpass, 
where brightness temperature is defined by equating the Planck function with the mean surface 
brightness. We note that our best-fit value for the depth of the transit (as characterized by ratio of 
the plantary to stellar radii) is slightly smaller than previous published values 6 ; this difference is 
most likely due to the effect of spots on the star. Large spots would increase the apparent depth of 
the transit at visible wavelengths, while having a minimal impact at 8 /xm. 
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Figure 3: Brightness estimates for 12 longitudinal strips on the surface of the planet. Data are 
shown as a colour map (a) and in graphical form (b); see below for details. We assume the planet 
is tidally locked, and we approximate it as being edge-on with no limb-darkening, so that the 
brightness of the 2th slice is Fi(sin0j )2 — sin0j i) where — tx/2 < ^,1,^,2 < tt/2 are the edges 
of the visible portion of each slice and Fi is the flux from a slice when it is closest to us. We bin 
the light curve into 32 bins with 4000 data points each, excising the data during the eclipses. We 
define our goodness-of-fit parameter as x 2 + A Y$L\ {F% ~ Fi-i) 2 > where x 2 is the goodness of fit for 
the light curve and the second term is a linear regularizing term which enforces small variations in 
adjacent slices for large A and allows a unique solution for Fj for a given value of A. We optimize 
this function using a 1000-step Markov Chain Monte Carlo method to determine the planetary 
flux profile and corresponding uncertainties. We chose a value for A that produced a reasonable 
compromise between the quality of the fit and the smoothness of the final brightness map. We 
varied both the size of the bins and the number of longitudinal slices and our resulting slice fluxes 
are robust. The brightness values in b are given as the ratio of the flux from an individual slice 
viewed face-on to the total flux of the star, with la errors, a is a Mollweide projection of this 
brightness distribution, with an additional sinusoidal dependence on latitude included (note that 
the data provide no latitude information). This plot uses a linear scale, with the brightest points in 
white and the darkest points in black. 
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